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We study the extended Anderson model on the two-dimensional Penrose lattice, combining the 
real-space dynamical mean-field theory with the non-crossing approximation. It is found that the 
Coulomb repulsion between localized and conduction electrons does not induce a valence transition, 
but the crossover between the Kondo and mixed valence states in contrast to the conventional 
periodic system. In the mixed-valence region close to the crossover, nontrivial valence distributions 
appear characteristic of the Penrose lattice, demonstrating that the mixed-valence state coexists 
with local Kondo states in certain sites. The electric reconstruction in the mixed valence region is 
also addressed. 


I. INTRODUCTION 

Quasicrystals have been receiving a lot of attention 
since its discovery [l|] . One of the interesting exam¬ 
ples is the rare-earth compounds Au-Al-Yb synthesized 
recently [2j. In the alloys, there are the quasicrys¬ 
tal AusiA^Ybis and the approximant Au 5 iAl 35 Ybi 4 , 
where the Yb icosahedra are arranged in quasi-periodic 
and periodic structures, respectively. The former qua¬ 
sicrystal shows non-Fermi liquid behavior with a non¬ 
trivial exponent in the specific heat and susceptibility, 
while the latter heavy fermion behavior m. This fact 
gives us a stage to study the role of the quasi-periodic 
structure in the strongly correlated electron systems 
7]. It has also been reported that the valence of the ytter¬ 
bium in these compounds is intermediate which indicates 
the hybridizations between 4/- and the conduction elec¬ 
trons [3]. This naturally arises a fundamental question 
how heavy fermion behavior is realized in the strongly 
correlated electron systems on the quasi-periodic lattice. 
An important point is that each lattice site in the quasi- 
periodic system is not equivalent, which is contrast to 
the conventional periodic system. Therefore, it is highly 
desired to clarify how the local geometry affects valence 
fluctuations and the Kondo state is realized by treating 
the quasi-periodic structure correctly. 

In this paper, we study valence fluctuations in the ex¬ 
tended version of the Anderson lattice model (EALM) [gj 
which includes the Coulomb interaction between the con¬ 
duction and /- electrons. We consider here the two- 
dimensional Penrose lattice (see Fig. [TJ as a simple quasi- 
periodic lattice. To take into account local electron cor¬ 
relations in the quasi-periodic lattice, we apply the real- 
space dynamical mean-field theory (RDMFT) mil to 
the model. We then study how valence fluctuations are 
affected by the quasi-periodic structure at finite temper¬ 
atures, comparing with the results for the EALM with 
the periodic lattice. The electric reconstruction in the 
mixed valence region is also addressed. 

The paper is organized as follows. In Sec. [Ill we in¬ 
troduce the EALM and briefly summarize our theoretical 
approach. In Sec. ED we study how valence fluctuations 



are enhanced in the system on the Penrose lattice. A 
brief summary is given in the last section. 


II. MODEL AND METHOD 

We study valence fluctuations in the EALM jg[ , which 
should be described by the following Hamiltonian, 

^ = (4 c i + h.c.) - V YJflcu, + h.c.) 
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where (i, j) denotes the summation over the nearest 
neighbor sites, Ci a ( fi a ) is an annihilation operator of 
a conducting electron (an /-electron) with spin <r(=t,i) 
on the ith site, n\ a = c\ a c i<7 and n{ a = f\ji a ■ Here, 
t is the hopping amplitude, V is the hybridization be- 
tween the conduction and / states, and e/ is the energy 
level of the / state. Uf ( U c f ) is the Coulomb interaction 
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in the / level (between the conduction electrons and /- 
electrons). 

The model with U c f = 0 has been studied by means 
of DMFT fin , where the competition between vari¬ 
ous phases has been discussed [13l - [l5l | . In the case of 
e/<0<6/ + [//, there should be one electron in the 
/-level, and the Kondo state is realized at low temper¬ 
atures. When the total number of electrons is close to 
half filling, the heavy-metallic Kondo state is realized 
with ( 77 /) ~ 1. In the case of €f ~ 0, the number of 
/-electrons is intermediate (( 77 /) < 1) and the mixed- 
valence state is realized. It is known that the crossover 
between the Kondo and mixed valence states appears 
by varying e/. The introduction of the interaction U c f 
enhances valence fluctuations and leads to stabilize the 
mixed valence state 0. At a certain critical point, the 
valence susceptibility diverges and the second order phase 
transition occurs between these two states (16]. Beyond 
the critical point, the valence transition is, in general, of 
first order, and a jump singularity appears in the number 
of /-electrons 77 /. This transition has been discussed in 
the EALM on the simple lattices such as one-dimensional 
chain (l6|, hypercubic 0 and Bethe lattices fl8j . 

Valence fluctuations in the quasi-periodic lat¬ 
tice, which should be important in the compound 
AusiA^Ybis, have not been discussed so far. To clarify 
this, we consider here the EALM on the two-dimensional 
Penrose lattice. In the lattice, a site is placed on each 
vertex of the rhombuses, as shown in FigHJ Since 
this lattice does not have a translational symmetry, 
lattice sites are distinct from each other. We note 
that the coordination number Z ranges from 3 to 7, 
except for edge sites. In infinite system ( N —>> 00 ), 
the average and its standard deviation are Z = 4 
and cr = ^TOr- 2 ~ 1.208 0- In the paper, we 
treat the Penrose lattice with the five-fold rotational 
symmetry, which is iteratively generated in terms of 
inflat ion-deflat ion rule [20|. In the noninteracting 
case Uf = U c f = 0, the Hamiltonian is numerically 
diagonalized. The densities of states (DOS) for the 
conduction and /- electrons are shown in Fig. (2J The 
tiny hybridization gap appears around u) = 0 and sharp 
peaks in the DOS for the /-electrons appear at the edges 
of the gap. In contrast, DOS of the conduction electrons 
are widely distributed, where the bandwidth W ~ 8.6L 
These are similar to those for the conventional periodic 
Anderson model. We also find the macroscopically 
degenerate states at uo = TV. These states are regarded 
as the bonding and antibonding states for the confined 
states discussed in the tight binding model on the 
Penrose lattice j2l|, [22|. These confined states may have 
little effects on low temperature properties in our model 
when we focus on the EALM with a fixed n = 1.9, where 

n = + n L)/ N - 

To discuss low temperature properties in the EALM 
on the Penrose lattice, we make use of the RDMFT 
[TH which takes local electron correlations into account. 
This treatment is formally exact for the homogeneous lat¬ 



FIG. 2. (Color online) DOS (upper panel) and integrated 
DOS (lower panel) for the conduction and / electrons in the 
noninteracting system with e/ = 0 and V/t = 0.5. 


tice in infinite dimensions, and enables us to obtain the 
reliable results if spatially extended correlations are neg¬ 
ligible. In fact, the method has successfully been applied 
to correlated systems such as surface 23|, interface a 
superlattice (25j| , ultracold atoms j26|, and topological 
insulating systems [27]. 

In the RDMFT method, the selfenergy should be site- 
diagonal as [S\ij c r = 'E, i(7 5ij, where Sij is a Kronecker’s 
delta function. The lattice Green function is then given 
by 



(iu n + /i) 1 



Uj 0 

0 0 


( 2 ) 

where 1 is the identity matrix, uj n = (2n + 1)ttT is the 
Matsubara frequency, fi is the chemical potential, and 
T is the temperature. The lattice model is mapped to 
effective impurity models dynamically connected to each 
heat bath. The effective imaginary-time action for the 
ith site is given as 



where = (c\ a fj a ) are Grassmann variables and 
Qq\t) is the Weiss effective field imposed on the self- 
consistency condition. The Weiss mean-field is obtained 
from the Dyson equation of the effective model, 

^(^n)- 1 = [Ga(^n)] 4 7 + S iCT . (4) 

When the RDMFT is applied to the EALM on the Pen¬ 
rose lattice, one solves the effective impurity models N 
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times by iteration. Here, we use the non-crossing ap¬ 
proximation mm as an impurity solver. This has an 
advantage in treating strong correlations at finite tem¬ 
peratures less expensive than the other numerical tech¬ 
niques such as the continuous-time quantum Monte Carlo 
method [32[ and numerical renormalization group [33j. 

To discuss how the valence are affected by electron cor¬ 
relations in the Penrose lattice, we calculate the number 
of /-electrons n{ a at the ith site. The valence suscepti¬ 
bility is defined by 


dnf 


Xv = ~ 


de f ’ 


(5) 


where nf = /N. In this paper, it is deduced in 

terms of the numerical derivative of the valence. In the 
following, we fix Uf/t = 80.0 and V/t = 0.5, and the 
total number of particles as n ~ 1.9. We set t — 1 as 
unit of energy . 


III. RESULTS 

In the section, we discuss how valence fluctuations de¬ 
velop in the correlated quasi-periodic systems. By per¬ 
forming RDMFT with the non-crossing approximation, 
we obtain the results for the system with U c f/t = 0,16 
and 36 at the temperature T/t = 0.2. We show in Fig. [3] 
the distribution of the number of /-electrons and valence 
susceptibility. When U c f = 0, the system is reduced to 
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FIG. 3. (Color online) Distribution of the number of /- 
electrons (a) and valence susceptibility (b) as functions of 
the energy of /-level in the system with N — 1591 at the 
temperature T/t = 0.2 when U c f/t = 0, 16 and 36. 


in the case U c f/t = 0, the quasi-periodic structure has 
little effect on low-temperature properties, at least, in 
the case with Uf/t = 80.0, V/t = 0.5 and T/t = 0.2. 

The introduction of the interaction U c f leads to the en¬ 
hancement of valence fluctuations in the system. In fact, 
the peak in the valence susceptibility becomes sharper 
when U c f increases, as shown in Fig. 03 Nevertheless, we 
could not find the divergence of the valence susceptibility 
even when the interaction is rather large (U c f/t = 36). 
The above results are contrast to that of in the EALM 
on the periodic lattice, where valence fluctuations are en¬ 
hanced around ef ~ ~U c f and the valence transition is, 
at last, induced. This discrepancy should originate from 
the geometry of the lattice. Namely, in the conventional 
periodic model, all lattice sites are equivalent. Therefore, 
the valence is suddenly changed at a certain point, where 
the first-order phase transition occurs between the Kondo 
and the mixed valence states. By contrast, in the Penrose 
lattice, each site is not equivalent, as shown in Fig. [U In 
our model, the bare onsite interactions Uf and U c f are 
uniform, but the coordination number Z depends on the 
site. Furthermore, the site geometry beyond the near¬ 
est neighbor sites is rather complex. Roughly speaking, 
local interactions are effectively modified, depending on 
the geometry around a certain site. This yields the site- 
dependent renormalization. Fig. [3] (a) shows the clear 
valence distribution in the crossover and mixed-valence 
regions. Therefore, we can say that the distribution of 
the effective interaction suppresses valence fluctuations 
and the crossover, instead, occurs between the Kondo 
and mixed valence states. 

By performing similar calculations, we determine the 
crossover line, where the valence susceptibility Xv has 
maximum. Fig. 0 ] shows the phase diagram for the sys¬ 
tem with Uf/t = 80 at the temperature T/t = 0.2. When 



the Anderson lattice model. In the case of ef/t < —3, 
the Kondo state is realized with 77 / ^ 1.0. As e/ ap¬ 
proaches 0, nf decreases and the mixed valence state is 
realized. Since no singularity appears in the curve, the 
crossover occurs between the Kondo and mixed-valence 
states. This is roughly determined by the maximum of 
the valence susceptibility. It is also found that the va¬ 
lence can be represented by a single curve although lat¬ 
tice sites are distinct from each other. This implies that 


FIG. 4. (Color online) Density plot of the standard deviation 
for the valence n/ in the system with T/t — 0.2 and N — 1591. 
The dashed line represents the crossover between the Kondo 
and mixed valence region. 


€f and U c f are small, the Kondo metallic state is real¬ 
ized with nf ~ 1. On the other hand, the mixed valence 
state is realized in the large U c f region. Furthermore, we 












4 


calculate the standard deviation of the valence to discuss 
site-dependent properties more clearly. The results are 
shown as the density plot in Fig. 0 ] In the Kondo regime, 
the valence is almost unity at each site, and thereby this 
quantity is negligible. On the other hand, in the mixed- 
valence regime (jif ^ 1), the value is finite and has the 
maximum around (e//t, U c f/t) ~ (—4,8), as shown in 
Fig. m 

To discuss how site-dependent properties emerge at fi¬ 
nite temperatures, we show in Fig. [5] (a) the temperature 
dependence of the valence in the system with €f/t = —5 
and Ucf/t = 10. At high temperatures T/t > 4, the sys¬ 
tem must be governed by the larger energy scales [//, U c f 
and €f. In the case, the valence little depends on the local 
geometry and a single peak appears in the cross-section of 
the valence distribution, as shown in Fig. [5](b). Decreas¬ 
ing temperature ( T/t = 1), a single peak is split into five 
peaks, as shown in Fig. [5] (c). This suggests that local 
geometry around a certain site affects low temperature 
properties, as discussed above. To proceed further discus¬ 
sions at lower temperatures, We employ de Bruijn’s nota¬ 
tion to classify the sites in the Penrose lattice [34| : eight 
kinds of vertices are denoted by D, J, Q, K , S3, S4, S and 
S5, as shown in Fig. [ 6 ] (a). Examining peak structures at 
T/t = 1 carefully, we find that the number of /-electrons 
nf ~ 0.22,0.19,0.17,0.15, and 0.14 are corresponding to 
the vertices S 3 , S 4 , {J, S, S 5 }, K , and {D,Q}, respec¬ 
tively. This means that at this temperature, the valence 
distribution depends on the coordination number rather 
than the vertex. Further decreasing temperature, many 
peaks appear in the valence distribution, and the vertex 
type becomes important. We find in Fig. [5] (a) that the 
valence for the sites with Z = 4 and 5 are divided into 
some classes at lower temperatures. As for the sites with 
Z = 5, the six kinds of peaks appear in the cross-section 
of valence distribution at T/t = 0.33, as shown in Fig. [5] 
(d). We find that these peaks originate from the local ge¬ 
ometry of the six vertices, as shown in Fig. [ 6 ] (b). Further 
decrease in temperature increases the valence for the ver¬ 
tices J-S5 and S-D, while it decreases for the vertices S5- 
J and J-S 3. Therefore, we can say that low temperature 
properties in the quasi-periodic system strongly depend 
on the local geometry. It is also found that each peak 
is divided into some small peaks at lower temperatures. 
This means that the longer range electron correlations 
become more important. Similar behavior appears for 
the sites with Z = 4. Namely, three kinds of peaks in the 
cross-section appears at lower temperatures, as shown in 
Fig. [5] (e). These originate from three kinds of vertices, 
as shown in Fig. [ 6 ] (c). Since there appears a difference 
in the configuration of the next-nearest-neighbor sites in 
three vertices, the small splitting appears even at lower 
temperatures, as shown in Fig. [5](e). 

To confirm that all vertices shown in Fig. [6] con¬ 
tribute the corresponding peaks, we independently count 
the number of vertices with Z = 5 [19| and Z = 4 in 
the Penrose lattice with N 00 , as shown in Tables U 
and HU The obtained probabilities are shown as the open 


TABLE I. The probabilities p of vertices with Z — 5 


Vertex S-D 

S5-J 

J-S3 

J-SA-a 

J-SA-/3 

J-S5 

p 1/a/5t 5 

1 /V$T 7 

1/r 6 

2/r 7 

1/r 7 

VZ/t 7 


TABLE II. The probabilities p of vertices with Z — 4 

Vertex S3-A-S3 S3-A-S4 S4-A-S4 
p 2/t 8 2/t 9 1/t 8 


circles in Figs. [5] (d) and (e). We find that these are 
consistent with the numerical results for the finite sys¬ 
tem N = 4181. This implies that the site-dependent 
renormalization characteristic of the Penrose lattice in¬ 
deed occurs. 

We would like to comment on the distribution of the 
valence. Fig. [5] (a) shows that the valences for the S3 
and S4 vertices rapidly increase, while for the D and 
Q vertices monotonically decrease as decreasing temper¬ 
ature. We note that the valence for the S3 vertices 
becomes much larger than the others and approaches 
unity, suggesting the formation of the local Kondo state. 
Therefore, we can say that at low temperatures, the 
mixed-valence sites with n{ < 1 coexists with the lo¬ 
cal Kondo sites with n{ ~ 1. To clarify how the above 
site-dependent properties affect the diffraction pattern, 
we calculate the quantity Ac = |rz^ | 2 , where 

( 6 ) 

i 

The results for the system with U c f/t = 10 at T/t = 0.2 
and 1.0 are shown in Fig. 0 When T/t — 1.0, the va¬ 
lence distribution is almost uniform and thereby A is 
reflected by the Penrose lattice. On the other hand, at 
the low temperature T/t — 0.2, an electric reconstruction 
is realized, where a rather large valence appears in the 
S3 vertices, as discussed above. This yields additional 
peak structures in /&, as shown in Fig. 0 Although this 
obtained pattern is characteristic of our model, the elec¬ 
tric reconstruction originating from electron correlations 
and local geometry is common to the quasi-periodic sys¬ 
tems. Therefore, we believe that such valence crossover 
should be observed experimentally in the quasicrystals 
with strongly correlated electrons. 

In the paper, we have studied the extended Ander¬ 
son model on the two-dimensional Penrose lattice. It 
has been clarified that the valence transition does not 
occur, but the crossover occurs with electric reconstruc¬ 
tions. Since we neglect various specific features for the 
compound Au 5 iAl 34 Ybi 5 such as the Tsai-type clusters 
and three-dimensional quasi-periodic structure, we could 
not explain the nature of the quantum critical behavior. 
It is an interesting problem to clarify which effect sta- 










5 




loos 




i>S3 J-S4-a| J-S4-0 J-S5 s_o 

A _ ."U ■ M. 


0.15 / 0.2 

rr 


Z=3 

Z=5 . 
Z=4 




0.15 r 05 055 


0.06 

(e) T<W05 — 

TJI-Q3S - . 

40.03 

P • 

A 

0 

mAA 


FIG. 5. (Color online) (a) The number of /-electrons as a function of the temperature T/t when Cf/t = —5, U c f/t = 10 and 
N = 4181. (b) and (c) show the cross-sections of the valence distribution at T/t = 10 and 1. (d) and (e) show the cross-sections 
at low temperatures with Z — 5 and Z — 4, respectively. 
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ometry in the Penrose lattice, which induces the electric 
reconstruction rather than the valence transition. In the 
mixed valence region close to the crossover region, the 
local geometry plays an important role in understand¬ 
ing low-temperature properties, where the mixed valence 
state coexists with local Kondo states. 
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FIG. 6. (Color online) (a) Classification of vertices in the Pen¬ 
rose lattice. The number represents the coordination number 
for each vertex. Detailed classification of vertices with Z — 5 
(b) and Z = A (c). 
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bilizes quantum critical phenomena in the quasicrystal, 
which is now under consideration. 


IV. CONCLUSIONS 

We have investigated valence fluctuations in the ex¬ 
tended Anderson model on the two-dimensional Penrose 
lattice, combining RDMFT with the non-crossing ap¬ 
proximation. Computing the number of /-electrons and 
the valence susceptibility, we have revealed that the va¬ 
lence transition does not occur in the vertex-type Penrose 
model even when the interaction between the conduction 
and /-electrons is large. This originates from the local ge- 
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FIG. 7. (Color online) Profiles of the quantity log 10 /k in the 
system (N — 4181) with €f/t — —5 and U c f/t = 10 when 
T/t — 0.2 (left panel) and T/t = 1.0 (right panel). 
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